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HYPERVELOCITY  IMPACT: 

DEPENDENCE  OF  CRATER  DIMENSIONS  ON  IMPACT  VELOCITY 

ABSTRACT 

Craters  in  copper  and  lead,  produced  by  hypervelocity  impact,  have 
been  measured  and  the  dimensions  correlated  with  impact  velocity.  The 
results  indicate  that  craters  scale  with  about  the  1.7  power  of  velocity, 
in  agreement  with  computer  physics  results  based  upon  hydrodynamic 
calculations. 
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INTRODUCTION 

i 

Numerous  empirical  correlations  of  crater  dimensions  with  projectile 
impact  velocity  have  been  presented  by  various  investigators.  In  general, 
these  correlations  can  be  classified  in  either  of  the  following  two  groups: 

1.  Crater  linear  dimensions  varied  with  the  one-third  power 
of  the  projectile  velocity  and  crater  volume  with  the  first  power  of 
velocity,  leading  to  a  statement  that  craters  scaled  with  the  projectile 
momentum. 

2.  Crater  linear  dimensions  varied  with  the  two-thirds  power 
of  the  projectile  velocity  and  crater  volume  with  the  square  of  the 
velocity,  leading  to  a  statement  that  craters  scaled  with  projectile 
energy. 

Actually,  it  has  long  been  recognized  that  both  points  of  view 
represent  oversimplifications  and  are  useful  only  for  purposes  of  inter¬ 
polation  and  modest  extrapolation.  It  was  at  first  assumed  that  the 
dependence  of  crater  dimensions  on  impact  velocity  could  not  be  arrived 
at  without  a  complete  treatment  of  the  crater  formation  process,  including 
the  later,  strength-dependent  portion.  This  supposition  is  too  restrictive, 
however.  Certainly,  the  size  of  a  given  crater  cannot  be  computed  by 
considering  only  the  hydrodynamic  portion  of  the  process,  without  treat¬ 
ing  strength  effects.  Yet,  certain  trends  can  be  discerned  Just  from 
performing  computations  in  the  regime  where  no  doubt  exists  as  to  the 

validity  of  the  hydrodynamic  assumption.  These  considerations  led  Walsh 
is 

and  co-workers  to  advance  the  concept  of  late-stage  hydrodynamic 
equivalence.  Simply  stated,  the  gross  features,  such  as  size  and  shape 
of  the  final  craters  from  two  different  impacts  will  be  the  same  if 
the  gross  transient  features  of  the  impacts  are  the  same  at  a  time  suf¬ 
ficiently  short  after  impact  so  that  pressures  in  the  target  are  large 
compared  to  the  strength  of  the  material.  That  is,  if  the  hydrodynamic 
solutions  of  the  flows  resulting  from  the  impacts  are  essentially  the 

-  ■ 

Superscript  numbers  denote  references  found  on  page  IS. 
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raw*  at  late  times  within  the  hydrodynamic  phase,  then  the  untreated 
post-hydrodynamic  phases  of  the  interactions  will  be  the  same.  This 
leads  to  a  scaling  of  craters  dependent  only  upon  the  mass  and  velocity 
of  the  projectiles  which  produced  them,  expressly  stated  by  the  relation: 


ml  V1 


3a 


“2 


3a 


Based  upon  consideration  of  momentum  and  density  distributions  in  the 
target  Walsh  and  co-workers  have  computed  3a  to  be  1.7^  -  0.06,  fro, 
both  particle-in-cell  and  continuous  Eulerian  formulations  of  two- 
dimensional  axisymmetric  problems. 


Observations  of  final  crater  dimensions  over  a  substantial  range  of 
projectile  velocity  can  be  used  to  determine  3a  experimentally.  This 
report  describes  this  determination,  summarizing  -erminal  crater 
observations  made  at  the  BRL  in  copper  and  lead. 


EXPERIMENTAL  OBSERVATIONS 

Chunky  projectiles  of  steel,  aluminum,  and  titanium,  were  accelerated 

from  several  types  of  explosive  devices  described  elsewhere.  These 

2  3 

include  the  BRL  Air-Cavity  Charge  ,  the  BRL  Inhibited— let  Charge  ,  and 

I4 

the  BRL  Micro-Particle  Accelerator  .  The  projectiles  struck  essentially 
semi-infinite  targets  of  lead  and  copper  at  velocities  up  to  12  km/sec., 
at  normal  incidence. 

Measurements  of  crater  dimensions  with  respect  to  the  level  of  the 
undisturbed  surface  were  made.  Crater  diameters,  normalized  by  dividing 
by  the  one-third  power  of  projectile  mass,  and  crater  volumes,  normalized 
by  dividing  by  projectile  mass,  are  plotted  versus  projectile  velocity 
in  Figures  1-U. 

These  data  were  then  fitted  to  functions  of  the  form: 

y  =  a  +  b  vc 
P 
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FIGURE  I  Norma li red  crater  diameter  data  for  normal  imnact* 
d?awP?hJouah8tJS'd%S  8  function  of  lBPact  velocity.  Th^line 
iap«t  5elo8it  h  da  COrresponds  t0  the  °-55  P°wer  of  the 


FIGURE  3  Normalized  crater  diameter  data  for  normal  impact 
in  lead  targets,  as  a  function  of  impact  velocity.  The  line 
drawn  through  the  data  corresponds  to  the  0.58  power  of  the 
impact  velocity. 


FIGURE  4  Normalized  crater  volume  data  for  normal  impacts  in 
lead  targets,  as  a  function  of  impact  velocity.  The  line  drawn 
through  the  data  corresponds  to  the  1.79  power  of  the  impact 
velocity.  K 


12 


by  applying  a  differential  correction  technique^  to  the  parameters  a,  b, 
and  c.  Thus,  in  order  to  determine  the  differential  corrections  da,  db, 
and  dc,  a  least  square  fit  was  performed  on  the  function 

dy  =  da  +  vc  db  +  b  vc  In  v  dc 
P  P  P 

The  differential  corrections  so  determined  were  then  applied  to  the 
original  estimates  of  a,  b,  and  c  and  the  process  repeated  as  many  times 
as  necessary  for  convergence  to  take  place.  Usually,  the  paramete  s 
converged,  to  five  significant  figures,  in  five  iterations. 

In  performing  the  fitting  operations,  velocities  have  been  ex¬ 
pressed  in  kilometers /second,  masses  in  gram,  diameters  in  centimeters 
and  volumes  in  cubic  centimeters. 


RESULTS 


The  results  of  the  analysis  described  above  are  summarized  in 
Table  I.  For  copper  targets,  drta  from  impacts  at  velocities  between 
2  and  12  km/sec  have  been  used.  For  lead  targets  data  between  0.05  and 
12  km/sec  have  been  used. 


Cu 


-0.165 


m 


3/3 


a 

0.311 


TABLE  I 
b  °b 

1.297  0 .2U7 


m 


— 1 . 1*3**  0.505  0. 1*61*  0.106 


co  a 

a  r 

0.5l*7  0.055  0.118 

1.767  0.091  1.532 


P,  D 

D  ~~/3  -°-0Q9  0.105  1.990  0.098  0.580  0.018  0.279 

m 

P 

V 

—  -0.238  0.376  1.796  0.097  1.791  0.025  U.5UU 

m 

P 

These  fits  are  plotted  in  Figures  1-1*. 
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For  the  diameter  fits,  c  represents  a  determination  of  a,  and  for 
the  volume  fits,  a  determination  of  3a.  Standard  deviations  are  shown 
for  all  parameters,  together  with  the  standard  deviation  of  the  residuals. 
The  intercept  a,  in  every  case,  is  not  sufficiently  far  from  the  origin 
to  be  significant. 


CONCLUSION 

The  experimental  values  of  a  (0.5^7  -  0.055  and  O.58O  -  0.018) 
and  3a  (1.767  -  0.091  and  1.791  -  0.025)  agree  quite  well  with  the 
computer  physics  results  of  0.58-  0.02  and  1.7^  -  0.06.  Thus,  it  can 
be  concluded  that  the  terminal  crater  data  do  not  contradict  the  hypothesis 
that  a  =  -  0.58  and  3a  =  -  1.7**. 
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